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拟 合 楼 层 反 应 谱 的 三 向 地 震动 时 程 生成 方法 研究 
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要 :为 获得 满足 规范 要 求 的 拟 合 楼 层 反 应 谱 的 三 向 地 震动 时 程 ， 六 阶 常 微 分 方程 的 本 征 函 
并 重 构 天 然 地 震动 时 程 ; 通 过 构建 影响 矩阵 ,和 迭代 计算 ee 数 的 幅 值 
系数 ,最 终 实 现 预 设 拟 合 精度 ;引入 正 交 化 方法 保证 三 0 并 通 
二 提供 验证 。 研 究 结果 表明 ,迭代 所 得 时 程 反应 谱 与 目标 楼 层 反 应 谱 的 拟 合 精度 较 高 , 迁 
过 程 中 时 程 反应 谱 均匀 一 致 地 向 目标 谱 逼 近 , 和 迭代 过 程 收 化 ,所 得 三 向 地 震动 时 程 均 无 震 球 现 
> 满足 规范 要 求 。 
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Abstract :In order to obtain the tri-dimensional seismic ground motion time histories which are compatible 
with the _ corresponding floor response spectra, the eigenfunctions of the sixth order ordinary differential 
equation are used as the basis functions to decompose and reconstruct the real recorded earthquake ground 
motion. By constructing the influence matrix, the amplitude coefficient of the eigenfunction is gradually 
adjusted through iterative calculation, and the preset fitting accuracy threshold is finally achieved. The 
orthogonalization procedure is introduced to ensure the statistical independence between the tri-directional 
time histories ,and two sets of numerical examples are provided. The examples show that it is capable of 


achieving high matching accuracy between the response spectra of generated time histories and the target 
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spectra. During the iteration ,the time history response spectrum uniformly approaches the target spectra. 


The iterative process converges; the generated seismic time histories are drift-free and satisfy the code 


requirements. 


Key words: seismic design and analysis; floor response spectrum; spectral compatible ground motion; 


cross correlation coefficient; eigenfunction 


核电 厂 的 抗震 分 析 中 ,往往 需要 3 个 正 交 方向 
(包括 2 个 正 交 水 平方 向 Hl 和 了 H2 ,和 1 个 竖 直方 向 
V) 的 地 震动 时 程 ,作为 激励 时 程 开展 动力 时 程 分 
析 ""。 由 于 核电 厂房 内 核电 设备 众多 且 构 筑 物 内 
部 结构 复杂 ,核电 设备 的 抗震 安全 性 分 析 成 为 核电 
厂 整 体 抗震 分 析 中 的 重点 和 难点 问题 “1 。 合 理 的 
地 震动 输入 是 核电 厂 内 部 重要 设备 (如 管道 ) 抗 震 
分 析 所 需 解 决 的 首要 问题 。 

楼 层 反 应 谱 是 核电 设备 抗震 输入 的 重要 依 
据 '" ,中国 住房 和 城乡 建设 部 发 布 的 《核电 厂 抗 震 
设计 标准 》( GB 50267 一 2019 ) ”中 国 国家 核 安 全 
局 《 核 设备 抗震 鉴定 试验 指南 》( HAF ' J0053 ) 5 、 
美国 核 管理 委员 会 NUREG- 0800《 标 准 审查 大 纲 》 
(SRP 3.7.1)5 加拿大 标准 协会 《核电 厂 抗震 鉴定 
分 析 指南 》( CSA N289. 3- 10) 等 相关 规范 均 指 
出 ,在 核电 厂 内 部 重要 构件 与 设备 的 抗震 动力 时 程 
分 析 中 ,所 输入 的 设计 地 震动 时 程 反应 谱 要 包 络 目 
标 楼 层 反 应 谱 , 并 满足 一 定 的 拟 合 精度 要 求 。 为 保 
证 三 向 地 震动 时 程 之 间 的 统计 独立 性 ,一 般 要 求 三 
向 地 震动 时 程 , 两 两 之 间 的 互相 关系 数 绝对 值 不 大 
于 0.165 0 或 0.39 1。 

当前 已 经 有 多 种 使 地 震动 时 程 反应 谱 与 目标 反 
应 谱 匹 配 的 方法 , 主要 可 分 为 基于 时 频 变 换 的 方 
法 ”在 时 域 琶 加 调整 函数 "中 以 及 调整 时 程 
分 量 幅 值 系数 的 方法 " ”1 。 然 而 ,由 于 楼 层 反应 谱 
在 形状 特征 上 与 普通 场地 反应 谱 相 比 有 一 定 特 殊 
性 2 , 现 有 技术 方法 中 尚 存在 以 下 不 足 : 个 拟 合 
精度 不 高 、 拟 合 过 程 有 “此 消 披 长 ”现象 ,由 于 时 域 
和 频 域 之 间 相 互 转换 的 复杂 性 ,以 及 传 里 叶 幅 值 谱 
与 反应 谱 之 间 的 非 线性 对 应 关系 ,使 得 基于 时 频 变 
换 的 方法 较 难 取 得 理想 结果 ;而 通过 在 时 域 局 部 鞋 
加 有 具有 特定 主 频 的 调整 函数 , 常 难 以 确保 不 对 其 余 
频率 点 的 反应 谱 值 产生 负面 影响 ,因此 , 拟 合 过 程 常 
不 收敛 ;@) 不 能 有 效 满 足 三 向 地 震动 分 量 间 的 统计 
独立 性 , 现 有 技术 少 有 仅 通过 一 次 计算 , 即 可 获得 满 
足 统计 独立 性 要 求 的 设计 地 震动 , 常 需 反 复 验 算 所 
得 时 程 。 


本 研究 基于 本 征 函 数 的 地 震动 时 程 分 解 ,构建 
影响 矩阵 反映 不 同 频率 地 震动 分 量 之 间 的 影响 关 
系 ,并 依 此 调整 不 同 频率 分 量 的 幅 值 系数 ,逐步 迭代 
以 使 时 程 反 应 谱 与 目标 楼 层 反 应 谱 的 拟 合 精度 满足 
要 求 。 通 过 在 迭代 流程 中 引入 正 交 化 方法 ,满足 不 
同方 向 分 量 之 间 的 统计 独立 性 要 求 。 本 研究 所 提出 
的 方法 计算 精度 高 ,迭代 收敛 性 好 ,同时 可 满足 规范 
中 相关 要 求 ,可 为 我 国 核 电厂 的 抗震 设计 与 分 析 提 
供 参考 。 


1 基于 本 征 函 数 的 地 震动 时 程 分 解 


为 获取 连续 的 加 速度 .速度 和 位 移 时 程 ,本 研究 
引入 六 阶 常 微分 方程 的 本 征 值 问题 ,通过 求解 得 到 
其 本 征 函 数 。 将 天 然 地 震动 时 程 利用 本 征 函 数 分 解 
并 重 构 ,可 得 到 与 原始 时 程 高 度 吻合 ,同时 可 避免 时 
程 零 味 现 象 的 一 组 加 速度 .速度 和 位 移 时 程 。 


1.1 六 阶 常 微分 方程 的 解 及 其 本 征 函 数 


考虑 如 下 关于 函数 p( 的 六 阶 常 微分 方程 的 
本 征 值 问题 
de hep =0,0<1<7 (1) 
满足 时 程 初始 与 终止 都 归 零 的 条 件 , 即 
p(0) = p(0) = (0) = p(7) = 
p(T) = p(T) =0 (2) 
可 得 到 其 通 解 为 
p(t) = Cicos(At) + Csin(At) + 


A At At 
e“ (Cscos + Casin 2 + 


2 2 (3) 


-a At .At 
A ee 
e (Cseos 7 + Cesin > ) 
对 上 式 分 别 求 一 阶 微分 和 二 阶 微分 可 得 
9(0) = on 交 


A A: .At 
Cse (V3cos > — sin | 十 
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“和 At . At 
学 
Gee (eos 7 + V3sin = 


C5e (V3eos SY + sin OY )- 


VE 


-3 {A 人 
Cee ( cos 了 + V3sin ”) (4) 
p(t) = -2Csin(At) -2C,cos(At) + 


Ba 
e? ee 3sin $4) + 


ci(V5cos 生 + sin 4)+ 


2 At . At 
2 Gs E 
© (eos 7 + V3sin 7 )+ 


eG (- V3cos 光 + sin 2) (5) 


2 
由 于 常 系数 Ci 中 > 有 非 零 解 ， 可 求 得 该 


问题 的 本 征 等 式 为 超越 方程 ， 人 的 根 或 本 征 
值 , 当 A，= (n+1)m/T, n=1, ,可 得 相应 


的 第 nn 个 本 征 函 数 gp, (1) 。 
图 1 为 持 时 30s 的 第 10 个 本 征 函 数 pio(b 及 
一 阶 和 二 阶 微分 示意 图 。 由 图 可 见 ,本 征 函 数 及 
其 一 、 二 阶 微分 的 中 间 形 态 均 与 三 角 函 数 非常 相似 ， 
具有 稳定 集中 的 频率 全 含量 ,上 且 两 端 始终 为 零 


0.6 


pid) 
So 


1! 


pid) 


PAD 4 


时 间 /s 
图 1 持 时 30s 的 第 10 个 本 征 函数 
Fig.1 The 10th eigenfunction with T=30s 


1.2 地 震动 时 程 的 本 征 函 数 分 解 


由 于 各 本 征 函 数 及 其 前 两 阶 微分 均 首 尾 为 堆 
因此 可 引入 本 征 函 数 为 基 函 数 ， 用 于 分 解 任意 持 时 
为 7, 且 首尾 归 零 的 天 然 地 震动 加 速度 时 程 4(1) 、 
速度 时 程 V(t) 和 位 移 时 程 P(t) ,1 € [0,7] 。 

首先 依据 地 震动 时 程 4(t1) 、V(t) 和 D(t) 所 涉 
及 的 频率 范围 [f,,,，f,,] ,用 下 式 求 得 本 征 频率 处 
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于 该 范围 内 的 本 征 函数 序号 为 
Na = {27 -11, No =L2Tf ~-1| (6) 
式 中 :运算 符号 [x| 表示 向 上 取 整 函数 ; Lx」 表示 
向 下 取 整 函数 。 
ee [f,,， fl 的 本 征 函 数 
Di (n= Nw， Naw +1,…,Ni) 作为 基 沪 数 , 按 
ee 


N 
Ninax 


A(t) = on p(t) (7) 


二 min 


利用 双 共 秋 梯 度 稳 定 法 ,可 求 得 构成 时 程 4(7) 
的 各 本 征 函 数 的 幅 值 系数 a ,n= N,N +1,…… 
Noo 基于 式 (3) ~ (5) ,考虑 加 速度 时 程 4(1) 、 速 
度 时 程 V(t) 和 位 移 时 程 p(t) 共用 相同 的 幅 值 系 
数 , 因 此 同时 可 得 


max 


Vt) = 5 08,(1) ,DG 


N 
N max 


= > api(t) (8) 
n= Nmin 


2 影响 矩阵 方法 


当 单 自由 度 受 非 平 稳 外 力 时 程 作用 时 ,一 般 可 
通过 卷 积 积分 求解 质点 的 全 响应 过 程 。 对 于 质点 质 
量 为 m, 自 振 圆 频率 为 o ,无 量 纲 阻 尼 系 数 为 《, 质 点 
系统 所 受 加 速度 时 程 激励 为 u(t) 时 ,质点 动力 方 
程 为 


X(t) + 2 X(t) + wx(t) = i,(t) (9) 
式 中 , x(t) 为 质点 相对 位 移 啊 应 时 程 。 当 阻尼 系数 
《“ 较 小 时 ,依据 卷 积 积分 原理 可 知 , 绝 对 加 速度 响应 


时 程 可 表示 为 
iD) ~ [insino(t -7)e dr 
二 (10) 
式 中 ,单位 肪 响应 函数 为 AD = er sea ， 
有 阻尼 圆 频 率 为 w = wV1 - 纪 。 


若 考虑 地 震动 激励 时 程 4(1) 利用 式 (7) 分 解 为 
本 征 函 数 的 表达 形式 ,通过 式 (10) 可 得 , 当 质 点 自 
振 圆 频率 分 别 为 ww， 普 = 1,2,… ,MM 时 的 加 速度 响 
应 时 程 亏 可 表示 为 


全 max 


Nn wa (hh, * A) 2 wn (h, > Cn 9,) 
五 = Nani 


min 


Ni ax 


二 > Qa 一 wh, * op,) (11) 


n= Nmin 
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式 中 ,单位 脉 串 响应 丽 数 为 ,= ee , 阴 尼 


圆 频 率 为 oo = wo V1 -4 。 那 么 , 单 自 由 度 体系 
加 速度 反应 谱 5,, ,可 根据 下 式 计算 


1 
ee ek 
i 六 max 


假定 加 速度 绝对 峰值 响应 发 生 于 时 刻 t = rw。 ， 
当 峰 值 响 应 (7) 三 0 时 , 取 é&，= 1 ; 当 峰 值 响应 
xn(Tn) < 0 时 , 取 é， =-1。 定 义 第 半 个 本 征 函 数 
ou = NanyNon+1 ,No ,在 频率 ww, 处 ,对 
于 反应 谱 值 5, 的 贡献 值 为 影响 因子 1 , 即 


Ls = ( 让 wh, 米 pa) | -， (13) 
那么 , 式 (12) 可 改写 为 
Na Nmax 
Sn 2 | py a — wih,, * 六 | = > Qn (én bn) 
ey max n= Nin 
(14) 


考虑 到 反应 谱 的 计算 共有 1 个 离散 频率 点 wn， 
= 1,2,…,M ,反应 谱 可 写 为 列 向 量 形式 $ = 
[SS Si ,那么 式 (14) 可 表示 为 如 下 的 矩阵 
形式 


S 二 la,a = [ ay, Qn max 

式 中 ,7 为 影响 矩阵 ,维度 为 Mx (N, -N +1)。 

假定 目标 楼 层 反 应 谱 $ 离散 于 M 个 频率 点 

9 M, 即 S = [SS ,Si Le 

那么 ,对 于 与 目标 谱 $ 相 匹配 的 地 震动 时 程 4(1) ， 

其 本 征 函 数 幅 值 向 量 4 可 通过 如 下 方程 求解 

S = 176 (16) 

由 于 所 用 反应 谱 离散 点 数 为 1 与 分 解 时 程 所 

用 本 征 函 数 数量 为 (NN, - Nu。+ 1) 往往 是 不 相等 

的 。 为 方便 求解 ,在 式 (16 ) 左 右 两 侧 乘 以 影响 矩阵 
的 转 置 7 , 即 


Qn ,MN 三 1 ， 2 


(T'Da = 1S (17) 

因此 ,可 求 得 线性 代数 方程 中 幅 值 向 量 a ,进而 

可 重 构 得 到 地 震动 加 速度 时 程 4(1) , 且 该 时 程 反应 
谱 比 原始 时 程 4(1) 的 反应 谱 更 接近 目标 谱 5 。 

需要 说 明 的 是 ,由 于 影响 矩阵 工 是 根据 峰值 响 

应 发 生 时 刻 7 计算 确定 的 , 当 目 标 谱 5 显著 区 别 于 

当前 时 程 反 应 谱 8 时 ,与 目标 谱 $ 匹配 的 时 程 4A(1) 

和 原始 时 程 4(1) 的 峰值 响应 发 生 时 刻 也 有 显著 差 
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异 。 因 此 ,本 研究 拟 基于 式 (17) 构 建 迭代 算法 ,使 
每 次 迭代 所 得 时 程 反 应 谱 逐 步 向 目标 谱 双 近 ,最 终 
实现 与 日 标 谱 的 匹配 精度 达到 任意 预 设 精度 。 


3 拟 合 楼 层 反应 谱 的 三 向 地 震动 时 程 


3.1 规范 对 于 拟 合 三 向 地 震动 时 程 的 要 求 


为 确保 抗震 设计 地 震动 时 程 的 合理 性 ,不 同 规 
范 中 ,针对 用 于 核电 厂 地 震 响应 分 析 的 地 震动 时 程 ， 
在 各 个 方面 均 给 出 了 明确 的 标准 。 通 过 整理 各 规范 
中 ,与 楼 层 反 应 谱 计 算 和 三 向 抗震 设计 地 震动 时 程 
相关 的 要 求 ,总结 了 几 项 广泛 适用 的 要 求 ,并 选用 较 
为 严格 的 参数 阔 值 ,作为 本 研究 中 地 震动 时 程 应 满 
足 的 要 求 ,以 期 通过 本 方法 计算 所 得 地 震动 时 程 能 
够 自然 适应 不 同 规范 ,具体 如 下 。 

1) 楼 层 反 应 谱 应 包含 两 个 水 平 正 交 方 向 的 谱 和 
一 个 竖 直 方向 的 谱 。 

2) 计算 楼 层 反 应 谱 时 ,频率 增 量 数值 应 按 规 范 
要 求 确定 。 

3 ) 用 于 结构 地 震 分 析 的 地 震波 , 需 为 真实 记录 
的 地 震波 或 由 真实 地 震 记录 调整 所 得 。 

4) 时 程 的 离散 时 间 间 隔 不 得 大 于 0.01 s。 

5) 设 计 地 震动 时 程 反应 谱 值 低 于 目标 反应 谱 的 
控制 点 数 不 得 多 于 5 个 ,在 所 有 控制 点 处 的 相对 误 
差 不 大 于 10% 。 

6) 每 条 设计 地 震动 时 程 的 强 震 持续 时 间 不 应 小 
于 6s, 总 持续 时 间 不 宜 小 于 25 s。 

7) 每 两 条 设计 地 震动 时 程 间 的 互相 关系 数 不 应 
大 于 0.16。 


3.2 拟 合 三 向 地 震动 时 程 的 步骤 


本 研究 以 本 征 函 数 分 解 地 震动 时 程 为 基础 , 引 
人 影响 矩阵 方法 构建 拟 合 楼 层 反应 谱 的 迭代 流程 ， 
并 引入 正 交 化 方法 确保 三 向 地 震动 时 程 间 的 统计 独 
立 性 ,具体 步骤 如 下 。 
1 ) 确定 水 平方 向 目标 楼 层 反应 谱 8 、S5 和 坚 
直方 向 目标 楼 层 反应 谱 5S 分 别 为 
Sa = [Sm Sm Shim], 
Sm = [Si,Sma，…,Stmor]， 
Sy = [SvSv…,Sv (18) 
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各 目标 谱 均 离散 于 W 个 频率 控制 点 处 凡 ， 
万 jy ,其 中 = fuin » fu = foro 

2 ) 选择 一 组 持续 时 间 为 7 了 的 种 子 地 震动 加 速度 
时 程 记 录 ( 包 含 正 交 的 水 平 Hl 和 H2 方向 分 量 ,和 
竖 直 V 方向 分 量 ) ,分 别 利用 式 (7) 和 (8 ) 将 其 分 解 
为 初始 选 代 时 程 4 (?) 、 A (5 和 4w (1) 。 


max 


A (1) = py 

AY (1) = en 

全 加 = 加 人 ) 2 0) (19) 
3) 设 置 反应 谱 拟 合 精度 。 对 于 目标 楼 层 反应 谱 


全 频率 范围 fj, <f < fi 内 ,满足 楼 层 反 应 谱 相 对 
误差 小 于 ei ; 对 于 结构 或 设备 重点 关注 频率 范围 
i 三 ffi 内 ,相对 误差 需 小 于 e, ;为 保证 工程 应 用 
中 重点 关注 频率 有 更 高 精度 要 求 ,通常 有 e，< el 。 

4) 求 解 水 平 Hl 方向 设计 地 震动 时 程 。 以 

A (+) 为 初始 迭代 时 程 ,利用 式 (17) 开 展 迭 代 计 
算 ， 直至 满足 预 设 的 拟 合 精 度 时 停止 。 可 得 与 目标 

谱 Sm 拟 合 的 Hl 方向 设计 地 震动 时 程 , 即 


Vmax 


4a(9 = DA (20) 


5) 求 解 水 平 H2 方向 设计 地 震动 时 程 。 以 
4 (1) 为 初始 迭代 时 程 ,在 第 i 次 迭代 中 
5a. 首先 按 下 式 利 用 Gram-Schmidt 正 交 化 调整 
水 平 卫 方向 时 程 与 水 平 Hl 方向 时 程 正 交 ( 即 互 相 
关系 数 为 0) , 则 


(4 44670) 
A -= AD (7) 加 HL A (1) 
HS H2 (A ,An) Hl 
Ninax 
> pa (21) 


式 中 ,符号 ”1 “表示 正 交 化 之 后 的 变 量 。 
5b. 利用 式 (17 ) 调整 地 震动 时 程 4 贞 (0b) 的 
幅 值 系数 使 反应 谱 向 目标 谱 副 近 , 得 到 AB (1) 。 
5c. 重复 步骤 5a 和 5b, 直到 满足 预 设 拟 合 精度 
为 止 。 可 得 与 目标 谱 Sm 拟 合 , 且 与 HIl 方向 正 交 的 
a 


Ain(t) = b, palt) (22) 


6) 求 解 竖 直 V 方向 设计 地 震动 时 程 . 以 
AW (0 为 初始 迭代 时 程 ,在 第 站 次 迁 代 中 
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6a. 按 下 式 利 用 Gram-Schmidt 正 交 化 调整 竖 直 

V 方向 时 程 ,使 其 与 水 平 Hl 和 H2 方向 时 程 的 正 交 
CAn 2) 

rn 

(A A 

i 


a 


-Fob (23) 


6b. 利用 式 (17) 调整 地 震动 时 程 Ay 0(t) 的 
幅 值 系数 使 其 反应 谱 向 目标 谱 副 近 , 得 到 AW (1) 。 

6c. 重复 步骤 6a 和 6b, 直 到 满足 预 设 拟 合 精度 
为 止 。 可 得 与 目标 谱 Sy 拟 合 , 且 与 水 平 Hl 和 H2 
方向 设计 地 震 时 程 正 交 的 竖 直 V 方向 设计 地 震动 时 
程 , 即 


N 
Ninax 


4v(t = 2 co pu 人 1) (24) 


7) 根 据 式 (20)、 (21) 和 (24) 所 得 加 束 束 度 时 程 ， 
利用 式 (8) 可 得 相应 速度 和 位 移 时 程 


N N 
N max N 


Vu (i) = > a.pa(t) ,Di (i) 2 anpn(t), 


全 三 


Nimax Ninax 
Viv(t) = > bip,(t) ,Din(t) = > b,p,(t), 
nS Won n= Nrnin 
Ninax Ninax 
Vt) = Dept) Dy) = cplt) 
5 nn 
(25) 


4 数值 算 例 


4.1 目标 楼 层 反应 谱 与 天 然 地 震动 


以 文献 [ 26 ] 提供 的 加 拿 大 某 核电 厂房 的 楼 层 
目标 谱 , 阻 尼 比 为 3% ,如 图 2 所 示 。 该 目 

谱 的 频率 下 限 与 上 限 分 别 为 上 =0.1 Hz 和 三。 
100 Hz。 

选取 1980 年 5 月 25 日 记录 于 Long Valley Dam 
(UprL Abut) 强 震 观 测 站 的 Mammoth Lakes 地 震动 
为 种 子 时 程 ,包括 两 个 正 交 的 水 平方 向 和 一 个 竖 直 
方向 的 地 震动 加 速度 、 速 度 和 位 移 时 程 记录 。 该 组 
地 震波 时 间 间 隔 均 为 0.005 s, 记 录 时 间 长 度 为 30s， 
每 一 条 时 程 有 6 001 个 数据 点 。Mammoth Lakes 地 
震动 为 正 斜 断层 ,震级 M, = 6. 06, 断 层 距 R= 
15. 46 km。 
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图 2 目标 楼 层 反 应 谱 (目标 谱 1) 
Fig.2 Target floor response spectra (targct spectra 1) 
利用 所 选用 的 三 向 种 子 地 震动 时 程 ,依据 本 研 
究 提出 的 方法 ,生成 分 别 与 三 向 目标 楼 层 反 应 谱 拟 
合 的 三 癌 设计 地 震动 时 程 。 为 验证 本 研究 方法 的 有 
效 性 与 精确 性 , 设 定 迭 代 时 程 反 应 谱 S(f) 与 目标 反 
应 谱 $ (六 的 相对 误差 按 下 式 计算 
_SCD -5'()) 
&( 力 三 gD x 100% 
依据 工程 结构 抗震 需求 设 定 如 下 精度 阀 值 。 
1) 在 全 频率 范围 内 , 拟 合 精度 e(j 不 大 于 5% 。 
2) 在 该 结构 重点 关注 频率 范围 [0.6,50]Hz 内 ， 
拟 合 精度 e(/ 不 大 于 1% 。 


(26) 


4.2 拟 合 结果 与 分 析 


考虑 目标 楼 层 反应 谱 的 频率 下 限 与 上 限 分 别 为 
fm =0.1 Hz 和 ,=100 Hz, 根 据 式 (6) 可 计算 得 到 
所 需 本 征 函 数 的 序号 范围 为 Wu =5 和 Vs。 = 
5 999 ,利用 所 选 本 征 函 数 重 构 天 然 地 震动 作为 初始 
迭代 时 程 。 

通过 影响 矩阵 方法 开展 迭代 计算 , 即 可 得 到 与 
楼 层 反 应 谱 高 精度 拟 合 的 三 向 地 震动 时 程 。 如 图 
3(a) 4(a) 和 5(a) 所 示 , 所 得 地 震动 时 程 反 应 谱 拟 
合 目 标 谱 的 精度 较 高 ,基本 可 完全 禾 盖 目标 谱 ; 且 从 
图 3(a) 图 4(a) 和 图 5(a) 中 灰色 的 中 间 迭 代 时 程 
反应 谱 可 见 , 随 着 迭代 的 进行 ,迭代 时 程 反 应 谱 均 匀 
一 臻 地 逐步 同 目 标 谱 荧 近 , 并 最 终 满 足 预 设 拟 合 精 
度 。 图 6 中 时 程 反应 谱 与 目标 谱 之 间 最 大 相对 误差 
随 着 迭代 计算 的 进行 逐步 减 小 ,反映 了 该 方法 的 有 
效 性 和 和 迭 代 收 敛 性 ,平均 每 次 迭代 用 时 约 7.0 s, 整 
体 计算 效率 较 高 。 从 图 中 速度 和 位 移 时 程 可 见 , 所 
得 地 震动 时 程 均 未 出 现 零 飘 现 象 , 且 各 向 地 震动 参 
数 均 符合 规范 要 求 ,具体 见 表 1。 
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图 3 水 平 HI 方向 反应 谱 拟 合 目 标 谱 1 的 
情况 和 所 得 设计 地 震动 


Fig.3 Compatibility between the response spectra of 


iterative time histories and the target spectra 1 along 


Fig. 4 


the Hl1 direction and the generated time history 
14 


1 水平 H2 方 向 目标 并 
wm 12[ 一 初始 迭代 时 程 
2 10| 一 中 间 时 程 
六 | 一 拟 合 时 程 
法 8 
i 
a 
: 4 
发 2 
0 
10 10. 10 10 
频率 /Hz 
村 (a) 地 震动 时 程 与 目标 谱 1 拟 合 情 况 
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图 4 水 平 卫 方向 反应 谱 拟 合 目标 谱 1 的 


情况 和 所 得 设计 地 震动 


Compatibility between the response spectra of iterative 


time histories and the target spectra 1 along the 


H2 direction and the generated time history 
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Fig.6 Relative error between response spectra of iterative 


time histories and the target spectra during iterations 
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表 1 拟 合 目标 谱 1 的 时 程 地 震动 参数 


Tab.1 Parameters of the generated time histories 


which are compatible with target spectra 1 


参数 水 平 Hl 方向 水 平 卫 方向 竖 直 V 方向 
euax/ 9 4. 828 4.963 4. 995 
enmin/ 92 —0.009 -0.016 —0.038 
emax_vip/ 9 0. 026 0.050 0. 005 
emin_vip/ 9 —0.009 -0.016 —0.003 
Tsup/s 9.675 10. 575 10.220 
Ni 4 4 3 
p pum =3.9x10- ,pl vy =1.3x107", 
Pi v =1.6 x10-1 


注 :eww enw 和 eww emiwiw 分 别 为 时 程 反应 谱 与 目标 楼 层 反 应 谱 的 
相对 误差 在 全 频率 [0.1,100]Hz 内 的 最 大 值 和 最 小 值 , 以 及 在 [0. 6， 
50] Hz 内 的 最 大 值 和 最 小 值 ; Towo 为 强 震 持续 时 间 , 取 Arias 烈度 自 
5% 至 75% 的 时 间 间隔 ;p 为 互相 关系 数 ,其 中 pvp 、 和 pi ,分 


别 为 三 向 地 震动 时 程 两 两 之 间 的 互相 关系 数 ;Ni 为 时 程 反 应 谱 连 续 
低 于 目标 谱 的 控制 频率 点 数 最 大 值 。 


5 案例 对 比 


5.1 目标 楼 层 反应 谱 与 天 然 地 震 即 


文献 [1 20] 提出 了 一 种 基于 HHT 的 多 维 地 震动 

拟 合 方法 ,本 节 将 以 该 文献 中 2 多 阻尼 的 三 向 楼 层 

应 谱 为 目标 (图 7) ,该 目标 谱 的 频率 下 限 与 上 限 
分 别 为 fi, =0.2 Hz 和 ,=100 Hz。 


40 
35 ---- 水 平 HI 方 向 ,-- 
各 了 2 
一 -一 坚 直 V 方 向 


反应 谱 /m 。 
GD 
So 


加 速度 
加 


10” 10° 10 10° 
频率 /Hz 


图 7 某 核 电厂 房 的 设计 楼 层 反应 谱 ( 目标 谱 2) 


Fig.7 Design floor response spectra for service building of 


a nuclear power plant (target spectra 2 ) 


采用 与 文献 [20] 相 同 的 地 震波 为 种 子 时 程 , 即 
在 Imperial Valley 地 震 (1979 年 ) 中 记录 于 Calexico 
Fire 台 站 的 三 向 地 震动 。 该 组 地 震波 时 间 间 隔 为 
0. 005 s, 记 录 时 间 长 度 为 37. 85 s, 每 一 条 时 程 有 
7571 个 数据 点 。Mammoth Lakes 地 震动 为 走 滑 断 
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度 设 置 与 4.1 节 相 同 。 


杨 兰 兰 ,等 : 拟 合 楼 层 反 应 
层 ,震级 ML, = 6. 53 ,断层 距 R,, = 10. 45 km, 拟 合 精 


5.2 拟 合 结 


， 


果 与 对 比分 析 


考虑 目标 楼 层 反应 谱 的 频率 下 限 与 上 限 分 别 为 


=0. 


n 


需 本 征 函 数 的 序号 范围 
7 。 利 用 所 选 本 征 函 数 重 构 天 然 地 震动 作为 初始 
迭代 时 程 。 利 用 本 研究 提出 的 方法 开展 欠 代 计算 ， 
分 别 经 过 710 .225 和 190 次 得 到 满足 精度 要 求 的 时 
程 。 时 程 反应 谱 拟 合 过 程 和 所 得 地 震动 时 程 如 图 8 
~10 所 示 , 所 得 设计 地 震动 时 程 满 足 规范 要 求 的 校 
核 结 果 如 表 2 所 示 。 


=100 Hz pe a 
为 N,, =15 和 WN 


2 Hz 和 六 


max 
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图 8 水 平 Hl 方向 反应 谱 拟 合 目标 谱 2 


的 情况 与 所 得 地 震动 时 程 


Fig.8 Compatibility between the response spectra of iterative 


time histories and the target spectra 2 along the 


HI1 direction and the generated time history 


从 所 得 时 程 反应 谱 与 目标 楼 层 反 应 谱 的 拟 合 精 


度 上 看 ,两 种 方法 调整 得 到 的 时 程 都 满足 规范 要 求 ， 
本 方法 所 得 时 程 反 应 谱 与 目标 谱 的 相对 误差 均 小 于 
1% , 相 较 于 文献 [ 20 ] 而 言 拟 合 精度 较 高 。 
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从 所 得 地 震动 时 程 波 形 上 看 ,本 人 研究 方法 所 得 
时 程 号 义 献 [20] 夺 入 导 时 程 均 无 球 移 现象 ,本 研究 所 


得 时 程 保留 了 更 多 天 然 地 震动 的 非 平稳 特性 。 


表 2 拟 合 目标 谱 2 的 时 程 地 震动 参数 


Tab.2 Parameters of the generated time histories 


which are compatible with target spectra 2 


参数 水 平 HIl 方向 水 平 H2 方向 坚 直 V 方向 
Cw/ 1.003 1.142 0.914 
emin/ % -0.019 -0.012 -0.003 
人 0.999 0.926 0.914 
emin_vip/ %% -0.019 -0.012 -0.003 
Tswp/s 10. 390 10.730 9. 995 
N, 3 5 5 
p Pip =1.8x10 ,pi =2.1x10™ 
pv =3.4x107 
40 yr pen 
--- 水 平 H2 方 向 目标 谱 
个 35| 一 初始 迁 代 时 程 
。 30 中 间 时 程 
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、,。| ~ 一 坚 直 V 方 向 目标 谱 
和 5[| 一 初始 迭代 时 程 
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Compatibility between the response spectra of iterative 


图 10 


Fig. 10 
time histories and the target spectra 2 along the 
V direction and the generated time history 
6 结 论 
核电 厂 及 其 内 部 设施 设备 的 抗震 性 能 是 核电 建 
设 和 发 展 所 关注 的 重要 内 容 。 为 获得 满足 规范 要 求 
的 拟 合 楼 层 反 应 谱 的 三 向 地 震动 时 程 作为 抗震 分 析 
的 动力 输入 ,本 研究 基于 本 征 函 数 分 解 天 然 地 震动 
时 程 ,提出 了 生成 拟 合 楼 层 反 应 谱 的 三 向 地 震动 时 
程 的 方法 ,并 以 加 拿 大 某 核电 三 楼 层 反应 谱 为 目标 
谱 ,提供 了 数值 算 例 。 所 得 主要 结论 如 下 。 
1) 本 研究 所 提出 的 方法 考虑 了 时 程 各 频率 分 量 
间 的 影响 ,以 及 各 频率 对 反应 谱 贡 献 的 影响 ,可 使 时 
程 反应 谱 与 目标 谱 全 频率 高 精度 拟 合 , 拟 合 过 程 均 
匀 且 收敛 。 
2) 本 方法 利用 本 征 函 数 首尾 自动 为 零 的 优点 ， 
使 得 本 征 函 数 分 解 后 的 地 震动 时 程 ,以 及 拟 合 上 日 标 
谱 的 地 震动 时 程 首尾 均 为 零 ,并 未 出 现 零 味 现 象 。 
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3) 本 方法 在 迭代 过 程 中 引入 正 交 化 方法 ,使 所 

得 拟 合 目 标 谱 的 三 向 地 震动 时 程 均 自动 严格 满足 互 

相关 系数 的 要 求 , 不 需 反 复试 算 或 重 算 。 
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